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Abstract 

Background: A fundamental pre-requisite for the clinical success in dental implant surgery is the fast and stable 
implant osseointegration. The press-fit phenomenon occurring at implant insertion induces biomechanical effects in 
the bone tissues, which ensure implant primary stability. In the field of dental surgery, the understanding of the key 
factors governing the osseointegration process still remains of utmost importance. A thorough analysis of the 
biomechanics of dental implantology requires a detailed knowledge of bone mechanical properties as well as an 
accurate definition of the jaw bone geometry. 

Methods: In this work, a CT image-based approach, combined with the Finite Element Method (FEM), has been 
used to investigate the effect of the drill size on the biomechanics of the dental implant technique. A very accurate 
model of the human mandible bone segment has been created by processing high resolution micro-CT image 
data. The press-fit phenomenon has been simulated by FE analyses for different common drill diameters 
(D A = 2.8 mm, D B = 3.3 mm, and D c = 3.8 mm) with depth L = 1 2 mm. A virtual implant model has been assumed 
with a cylindrical geometry having height L = 1 1 mm and diameter D = 4 mm. 

Results: The maximum stresses calculated for drill diameters D A , D B and D c have been 12.31 GPa, 7.74 GPa and 
4.52 GPa, respectively. High strain values have been measured in the cortical area for the models of diameters D A 
and D B/ while a uniform distribution has been observed for the model of diameter D c . The maximum logarithmic 
strains, calculated in nonlinear analyses, have been e = 2.46, 0.51 and 0.49 for the three models, respectively. 

Conclusions: This study introduces a very powerful, accurate and non-destructive methodology for investigating 
the effect of the drill size on the biomechanics of the dental implant technique. 

Further studies could aim at understanding how different drill shapes can determine the optimal press-fit condition 
with an equally distributed preload on both the cortical and trabecular structure around the implant. 



Background 

The use of dental implants has rapidly evolved since the 
advent of osseointegration, progressively replacing 
removable dentures in treatment of partially or com- 
pletely edentulous patients. A fundamental pre-requisite 
for the clinical success in dental implant surgery is the 
fast and stable implant osseointegration. 

The press-fit technique, which is adopted to ensure a 
primary stability condition for endosseous implants, 
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requires the diameter of the hole drilled within the jaw 
bone to be smaller than the implant major diameter. 
The bone tissues are affected by a biomechanical 
phenomenon, which is characterised by the mechanical 
properties of bone, implant materials, difference between 
implant and hole diameters and surrounding bone 
morphology. 

Recently, image-based approaches combined with 
Finite Element Analyses (FEA) have allowed effective 
stress-strain investigations in dental implantology. Den- 
tal implants can be virtually positioned within realistic 
models of human jaws reproduced from high definition 
CT image data with respect of the anatomical- 
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Figure 1 Micro-CT scan data. Micro-CT scan data of two different 

sections of the mandibular bone segment (35 urn resolution). 
L J 



physiological structures of bones. Worldwide, scientists 
have focused on this topic, especially to improve the suc- 
cess of endosseous implants [1-4]. A common goal is to 
understand the key factors of osseointegration processes 
following implant surgeries. Some researchers have 
investigated micro-displacements occurring at the bone- 
implant interface, while others studies have considered 
the load transfer at the interface to be more important 
in determining the correct mechanical stimulation of the 
osteoblasts, which are assumed to be responsible for 
bone tissue regeneration and the consequent osseointe- 
gration of the implant [5,6]. Generally, trabecular micro- 
structures of bones are modelled as homogeneous 
entities with particular mechanical properties and con- 
tiguity assumed at implant-bone interfaces. The contigu- 
ity conditions do not allow relative motions between the 
parts generating a continuum of stress distribution at 
the interface, where stresses are generally concentrated. 
Limbert [7] has considered the trabecular microstructure 
of the mandible bone and the discontinuity at the im- 
plant-bone interface by a Finite Element Analysis. Fur- 
ther studies investigate the preload condition generated 
by the insertion of the abutment screw in the implant 
for different designs of the screw-abutment system [6]. 




Figure 2 Edentulous bone segment. Different views of the 
edentulous bone segment generated in ScanlP. 




Figure 3 Drill CAD model. Positioning process of the drill CAD 
model. 

v_ J 



Nowadays, the selection of drill diameter with respect to 
the implant geometrical configuration is still done with- 
out any scientifically tested criteria. Natali [8] has ana- 
lysed the press-fit phenomenon occurring in oral 
implantology using the FE approach. In this study, the 
mandibular bone has been reconstructed by using CT 
data and attributing different mechanical properties in 
the cortical and trabecular regions. 

However, a more detailed study should also consider 
the trabecular microstructure of real bone tissue and the 
contacts associated to the relative movement between 
implant and bone. 

In this paper, an accurate model of human mandible 
bone segment is created processing high resolution 
micro-CT data by using image-based tools. The bio- 
mechanics of press-fit phenomena has been analysed by 
FE methods for different drill diameters. 
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Methods 

In this paper, edentulous bone segments of the right molar 
mandibular area have been acquired using a high resolution 
microtomography machine, the SkyScan micro-CT 
(SKYSCAN, Kartuizersweg, Kontich, Belgium). The sliced 
images have been generated with an isotropic resolution 




BEGINNING OF STEP 



MIDDLE OF STEP 



END OF STEP 



Figure 5 Interference fit algorithm. Contact constraint imposed at 
the bone-implant interface. 



of 35 Lim. Figure 1 shows the computed tomography 
images of two bone slices. 

The images from the micro CT have been processed to 
obtain the 3D models of the bone segments using the soft- 
ware ScanIP (Simpleware Ltd., Innovation Centre, Rennes 
Drive, Exeter, United Kingdom). 2D slice images have 
been processed by a methodology articulated in two main 
steps: segmentation and voxel creation. Initially, data have 
been converted to grey scale images. The region of interest 
has been isolated by thresholding the images and flood 
filling the masked area to avoid disconnected zones. The 
flood filling algorithm is based on a region-growing algo- 
rithm also known as paint bucket The obtained segmen- 
ted slices have been combined to produce isotropic 
voxels. The generated model is shown in Figure 2. The 
use of high resolution micro-CT data has allowed a 
detailed description of the trabecular bone microstruc- 
tures with a voxel size of 130 |im. 

The drill has been modelled assuming a cylindrical 
shape with different diameter values: D A = 2.8 mm, D B = 
3.3 mm, and D c = 3.8 mm and length L = 12 mm. The 
virtual implant model has been assumed with a cylin- 
drical geometry having height L = ll mm and diameter 
D = 4 mm. The drill models have been positioned into 
the bone model throughout the use of ScanCAD (Sim- 
pleware Ltd., Innovation Centre, Rennes Drive, Exeter, 
United Kingdom), which allows the placement to be 
driven by the corresponding 3D grey-scale micro-CT 
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Figure 6 Von Mises stress distributions. Top view of the Von 
Mises stresses induced by bone tissue expansion for models A, B 
and C, respectively. 



proposed in literature [5,10]. The implant, typically 
made of Titanium (E = 115 GPa), has been modelled as a 
rigid body considering the huge difference between the 
bone and implant Young s modulus and the effective re- 
duction of the simulation computational cost. 

The bone behaviour induced by inserting an implant 
with a major diameter into a smaller hole has been 
simulated using a contact algorithm (Interference Fit) 
based on the relative geometries, the material properties 
and the contact parameters. The assembly created after 
the implant positioning generates a volumetric overlap- 
ping between implant and bone shapes. 

The contact constraint, imposed at each constraint lo- 
cation (bone-implant interface) to avoid overlapping 
zones, can be considered as the current penetration (h). 
Penetration exists when h is positive (Figure 5). This 
constraint can be updated by specifying an allowable 
interference (v), which is ramped down over the course 
of a step. The specified allowable interference modifies 
the contact constraint as h - v = 0. 

The contact interaction property has been defined by 
a contact in the normal direction of the interaction (no 
penetration is allowed between the surfaces) and a fric- 
tion coefficient u = 0.2 in the tangential direction. 

The lateral sides of the bone segment have been con- 
strained in order to simulate symmetry conditions. The 
simulations have been run using Comsol Multiphysics 
Analysis and the Static step option. The static analysis 
approach has been used to analyse the stresses generated 
by the bone tissue stabilization following the press-fit 
phenomenon. 

In accordance with the Declaration of Helsinki (Hel- 
sinki, 1964), written informed consent was obtained 
from the patient and the study protocol was approved 
by the local Ethics Committee. (University of Sassari 
no. 985) 



images and taking into account the bone structure 
(Figure 3). The result of the virtual drilling process is 
shown in Figure 4. 

The edentulous drilled bone has been automatically 
meshed using the software ScanFE (Simpleware Ltd., 
Innovation Centre, Rennes Drive, Exeter, United King- 
dom). The meshing process has been carried out using 
EVoMaCs approach [9], providing meshes of approxi- 
mately 3 million elements (eight-node brick elements 
and four-node tetrahedral elements) and 900,000 nodes. 
The numerical analyses have been carried out by using 
the Comsol Multiphysics Analysis simulation software. 
The high 3D model accuracy has provided meshes free 
of any relevant volume/topology data loss. 

The bone tissue has been modelled using an isotropic 
linear elastic material (Young s modulus E = 15 GPa, 
Poissons ratio v = 0.3 and density p = 731 mg/cm 3 ) as 



Results and discussion 

The finite element analysis has been conducted for the 
three different drill diameter values. Figure 6 show the 
Von Mises stress distributions from a top view for mod- 
els A, B and C, respectively. In Figures 7 and 8, the Von 
Mises stresses are presented along the bucco-lingual and 
mesio-distal planes. In models A and B, the press-fit 
phenomenon is almost completely supported by the cor- 
tical structure, though the implant is in contact with both 
the cortical and the trabecular bone. The stresses on the 
trabecular zone are minimal compared to those charac- 
terising the cortical area. In model C the stresses are 
equally distributed across the cortical and trabecular 
structure reducing the generation of critical and localised 
loads. A uniform stress distribution, as reported in model 
C, is the optimal condition for the osseointegration 
process allowing a distribution of bone cell stimulation 
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Buccal Side < ► Lingual Side Lingual Side < ► Buccal Side 




Figure 7 Von Mises stresses distribution. View of the Von Mises stresses distribution on opposite sides of the bucco-lingual virtual cut. 



all around the implant. High stress concentrations could 
determine bone cell necrosis or lack of homogeneous 
osteoblast stimulation, which are both possible causes of 
implant destabilisation. The maximum stresses in models 
A, B and C are 12.31 GPa, 7.74 GPa and 4.52 GPa, 
respectively. 



According to the stress distribution, high strain values 
are measured in the cortical area for the models A and 
B, and a uniform distribution is observed in model C. 
The maximum logarithmic strains, calculated in non- 
linear analyses, are e = 2.46, 0.51 and 0.49 for models A, 
B and C, respectively. The logarithmic strain 
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distributions are presented in Figure 9 along the bucco- 
lingual and mesio-distal planes. 

The press-fit phenomenon indicates a mechanical con- 
nection of two structures based on contact pressure. 
However, a more appropriate definition, "interference- 



fit", could be used considering that each component mu- 
tually interferes to generate the deformation. 

Generally, the actual viscoelastic behaviour of the bone 
tissue imply that the elastic recoil of bones decreases over 
time. The bone deformation and the elastic recoil depend 



Frisardi et a I. Head & Face Medicine 201 2, 8:1 8 
http://www.head-face-med.eom/content/8/1 /1 8 



Page 7 of 9 




Frisardi et a I. Head & Face Medicine 201 2, 8:1 8 
http://www.head-face-med.eom/content/8/1 /1 8 



Page 8 of 9 



on the density of bone tissue, and the cortical bone has a 
less viscoelastic behaviour than cancellous bone [11,12]. 
The results obtained in the present work demonstrate that 
the cortical bone is stressed and deformed significantly, 
even with moderate values of interference -fit. Figures 6,7 
and 8 show higher von Mises values in the juxtacortical 
areas than in the bone marrow. 

A viscoelastic behaviour also depends on the degree of 
initial elastic deformation of the bone when the press-fit 
is achieved. The "assembly strain" (i.e., the deformation 
of the bone tissue upon the introduction of the implant) 
has an effect on both the contact pressure between bone 
and implants (e.g. the initial press-fit) and the rate at 
which the press-fit is dissipated through deformation 
and elastic relaxation [13,14]. 

The results shown in Figures 7, 8, 9 indicate that the 
only press -fit value (pre-load) evenly distributing stresses 
caused by the bone- implant contact is generated with 
model C; whereas models A and B yield a uneven stress 
distributions. However, the stresses occurring in the 
model C may not yet be high enough to yield the neces- 
sary press-fit persisting in order to assure the primary 
stability. Once the initial press -fit is dissipated, the im- 
plant could move in the bone under load. This could 
suggest the use of an "undersized" approach that would 
only accentuate an elevated cortical bone-implant load 
using the techniques at our disposal. 

The period between the initial press-fit and the healing 
or bone integration (intermediate period) is typically 
characterised by a relaxation of press-fit, which could be 
so fast to cause implant instability. The bone ingrowth, 
which leads to bone integration, is a very susceptible 
biophysical process, especially in the case of immediate 
loading of full arches where the absence of periodontal 
pressure receptors and the changes in the modulation of 
the trigeminal reflexes are no longer efficient. 

The results of this study suggest the development of 
an alternative surgical procedure able to determine uni- 
form stress and strain values between cortical and can- 
cellous bone. This method should not be classified as a 
generic "interference-fit", but as "differentiated interfer- 
ence-fit". 

A recent study [15] demonstrated that the undersized 
drilling technique, in addition to enhancing primary im- 
plant stability, might also achieve a translocation of bone 
particles. This biological phenomenon, comparable to 
the autologous bone graft, has a positive influence on 
the osteogenic response. However, this study shows that 
the difference between the implant diameter and the im- 
plant drill diameter is ~ 0.6 mm corresponding approxi- 
mately to model B. 

According to the most recent studies focusing on pri- 
mary stability [16] with an undersized approach [17], 
which is still equivalent to a press -fit of model B of this 



study, the occlusal load should also be considered, be- 
sides the obvious biophysical properties of the cortical 
bone and marrow. This does not correspond to a con- 
stant static and axial force, but to a dynamic, angular 
moment and extra-axial force with six degrees of 
freedom. 

Conclusions 

This study introduces a powerful, accurate and non-de- 
structive methodology for investigating the effect of the 
drill size on the biomechanics of the dental implant tech- 
nique. Further studies could aim at understanding how 
different drill shapes can determine the optimal press-fit 
condition with an equally distributed pre-load on both the 
cortical and trabecular structure around the implant. 
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